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We report experimental and analytical work on spin-transfer torque induced vortex dynamics 
in metallic nanopillars with in-plane magnetized layers. We study nanopillars with a diameter of 
150 nm, containing two Fe layers with a thickness of 15 nm and 30 nm respectively, separated by a 
6 nm Ag spacer. The sample geometry is such that it allows for the formation of magnetic vortices 
in the Fe disks. As confirmed by micromagnetic simulations, we are able to prepare states where one 
magnetic layer is homogeneously magnetized while the other contains a vortex. We experimentally 
show that in this configuration spin-transfer torque can excite vortex dynamics and analyze their 
dependence on a magnetic field applied in the sample plane. The center of gyration is continuosly 
dislocated from the disk center, and the potential changes its shape with field strength. The latter is 
reflected in the field dependence of the excitation frequency. In the second part we propose a novel 
mechanism for the excitation of the gyrotropic mode in nanopillars with a perfectly homogeneously 
magnetized in-plane polarizing layer. We analytically show that in this configuration the vortex can 
absorb energy from the spin-polarized electric current if the angular spin-transfer efficiency function 
is asymmetric. This effect is supported by micromagnetic simulations. 



I. INTRODUCTION 

Since the theoretical prediction [TJ [5] and experimen- 
tal demonstration 3, 4 J of spin-transfer torque induced 
effects many studies focussed on phenomena such as 
current induced magnetization dynamics and switching. 
Spin-transfer devices are promising candidates for future 
information technology. Current-induced switching be- 
tween discrete magnetization states may be employed in 
non- volatile memory such as the magnetic random access 
memory. As spin-polarized currents can propel steady 
spin precession [5;, spin-transfer torque devices are also 
envisaged to be used as integrated microwave sources. 
Therefore, finding highly tunable spin-transfer nano os- 
cillators (STNOs) is a matter of great current interest. 
After the discovery that spin-transfer torque can also 
drive the oscillatory motion of a magnetic vortex [6] , the 
gyrotropic vortex motion is considered as a promising 
STNO candidate. The magnetic vortex structure ap- 
pears at certain dimensions as the ground state of fer- 
romagnetic disks. The magnetization circulates around 
a core with magnetization perpendicular to the plane of 
the disk. The orientation of the magnetization in the 
core is called the polarization and the sense of rotation of 
the in-plane circulation constitutes the vortex chirality. 
This configuration is a consequence of the competition 
between the dipolar and exchange energies. Various ma- 
terial combinations have been used in previous studies 
on spin-transfer torque. In case of metallic spin valves, 
for example, Co/Cu systems are very common. In our 
study we use a molecular beam epitaxy-grown system 
containing Fe as ferromagnetic and Ag as spacer mate- 
rial [Fe/Ag(001)]. It has been calculated [7] that the 
Fe/Ag(001) interface exhibits a resistance that strongly 



depends on the spin direction thus providing for a high 
spin polarization. Additionally it has been shown [5] that 
an asymmetry in (spin channel averaged) ferromagnetic 
element and spacer resistances leads to an angular asym- 
metry in giant magnetoresistance (GMR) and also in spin 
torque efficiency. Experimental evidence for the appear- 
ance of this feature in Fe/Ag/Fe(001) spin valves has 
been provided by Lehndorff et al. [9] . Furthermore, it has 
been shown for the case of a thin (2 nm) extended polar- 
izer that vortex motion can be excited in Fe/Ag/Fe(001) 
pillars [TU] , 

In this work both magnetic layers have a finite lateral 
extension, which makes micromagnetic simulations of the 
system more feasible. As will be shown the measurements 
agree very well with the calculations. After the experi- 
mental part we will address the question how the vortex 
motion in a pillar with in-plane magnetized polarizer is 
excited and propose a novel mechanism. We show by an- 
alytical calculation and micromagnetic simulations that 
the angular asymmetry in spin-transfer torque efficiency 
can lead to energy absorption of the magnetic vortex from 
the spin-polarized electric current. 



II. SAMPLE FABRICATION 

After cleaning and annealing the GaAs(lOO) substrate, 
the sequence of metallic layers Fe (lnm) / Ag(150nm) 
/ Fe(30nm) / Ag(6nm) / Fe(15nm) / Au(50nm) is de- 
posited by molecular beam epitaxy. The 1 nm Fe layer 
serves as a seed layer for the adjacent Ag buffer layer, 
which will later on be the bottom electrode. It follows 
the spin valve layer system with an Au capping layer on 
top, which will also provide a part of the top electrode. 
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FIG. 1: Resistance vs. field dependence for low bias current 
(1 mA). The graphs are vertically shifted for better visibility. 
Top (bottom) graph: Applied magnetic field aligned with easy 
axis 1 (2). Insets: Schematic top views of the sample, the ar- 
row indicating the positive field direction. The measurements 
show that the two axes are not equivalent. 



All layers grow in a singlecrystalline manner [IT]. Both 
Fe layers have bcc structure and exhibit cubic magne- 
tocrystalline anisotropy. Next the areas of the prospec- 
tive bottom electrodes are covered with resist by optical 
lithography. The uncovered material is then removed 
with ion beam etching (IBE). Afterwards resist dots of 
150 nm diameter are placed into the 15 x 15 /xm 2 contact 
area located in the middle of the previously defined struc- 
tures by means of electron beam lithography. We employ 
a HSQ electron beam resist (fOx-12). Successively, the 
layers are milled down practically onto the bottom elec- 
trode by IBE resulting in pillars with diameters of about 
150 nm. We again use HSQ to insulate the nanopillars 
and planarize the sample. Additional protection and 
insulation is provided for by a 50 nm thick SiaN4 layer 
which is deposited by means of plasma-enhanced chemi- 
cal vapor deposition (PECVD). In the next step the sam- 
ple is covered with optical resist apart from 10 x 10 /xm 2 
contact windows above the nanopillars, which are opened 
by IBE so that the pillar tops are cleared from insulating 
material. A short reactive ion etching (RIE) step is addi- 
tionally performed to selectively remove remaining HSQ. 
Finally, the top electrode is deposited by a negative op- 
tical lithography step and subsequent lift-off. 
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FIG. 2: (a) Resistance vs. field dependence for a large bias 
current of —17 mA (electron flow from bottom to top elec- 
trode). The field is aligned with easy axis 2. Stars indicate 
high frequency excitations, (b) Spectra taken at each field 
value of the field sweep from positive to negative values (red 
curve in Fig. [2^a). (c) Enlargement of the vortex dynamics 
occurring between 197 and 86 mT. 



III. VORTEX DYNAMICS: EXPERIMENT AND 
SIMULATION 

All measurements are carried out at room tempera- 
ture. To characterize the dc behavior of the sample 
we apply a constant current while sweeping the mag- 
netic field and measuring the resistance (2-point measure- 
ments) . For measurements of the RF voltage we split the 



ac and the dc part of the voltage by a bias-T (frequency 
range: 45 MHz-26.5 GHz). The signal is then amplified 
by a 30 dB amplifier (0.5-26.5 GHz) and fed into a 50 GHz 
spectrum analyzer. The sample holder is home-made and 
contacts to the sample leads via spring pins with a diam- 
eter of 0.3 mm. The angles between the in-plane easy 
axes of the cubic anisotropy and the long axes of the top 
and bottom leads are 45° (insets in Fig. [l}. (The two in- 



3 







B 






L 







U i ' I ■ I ' — — ■ I 1 ' 1 

-600 -400 -200 200 400 600 



FIG. 3: Resistance vs. field dependence obtained from mi- 
cromagnetic simulations. The field is aligned with easy axis 
2. The symbols display the associated magnetization states. 
The simulation takes into account the Oersted field, which is 
derived from an inhomogeneous current density of —14 mA 
total current. The current distribution reflects the pillar and 
lead geometry. 



plane hard axes of the cubic anisotropy are parallel to the 
top and bottom leads respectively.) Fig. [I] shows mea- 
surements of the dc resistance at low currents for both 
easy axes (1st axis: black, 2nd axis: red). The difference 
of the two cases is likely to be caused by sample imperfec- 
tions (e.g. electrode contacting, deviations from perfect 
cylindrical shape of the nanopillar) which can lead to 
nonsymmetric current distributions and additional con- 
tributions to anisotropy. At small applied fields, most in- 
terestingly, the axis-2 measurement exhibits states with 
resistances as low as the saturation resistance. Instead 
we would expect intermediate up to maximum resistance 
values due to the dipolar antiferromagnetic coupling be- 
tween the layers. To investigate the sample characteris- 
tics under high current densities we apply a current of 
— 17 mA (Fig. |2^) (electron flow from bottom to top of 
the pillar) . We expect increasing influence of spin torque 
effects and the Oersted field on the disk magnetizations. 
The field is applied parallel to easy axis 2 and swept from 
negative to positive (sweep 1, black) and back to negative 
values (sweep 2, red). Both sweeps show an increase of 
the GMR as the applied field is driven from the satura- 
tion range towards zero. Once a sufficiently low external 
field is reached (— 51mT for sweep 1, 86 mT for sweep 
2) the sample enters a state of low resistance which per- 
sists up to large opposite field values (274 mT for sweep 
1), where for both sweeps we see abrupt steps in resis- 
tance. Afterwards the resistance decreases as the sample 
is driven into saturation. 

Micromagnetic simulations have been performed with 
the finite element code called TetraMag [12] , The mag- 
netization configurations in each disk for the different 
applied field values are provided by the micromagnetic 
simulations. Thus a clear physical understanding of 



the GMR signal is possible. Typical material parame- 
ters of iron, = 2.15 T (saturation magnetization), 
exchange constant A = 2.1 x 10~ n J/m and anisotropy 
constant K c = 48kJ/m 3 are used. The sample volume 
is discretized into irregular tetrahedrons with cell size of 
about 1 nm. The magnetic field is applied along the 2nd 
easy axis and linearly increased/decreased over a 40 ns 
simulation time between the two extrem values. The 
damping factor is set to a = 0.1. 

The simulations take into account the Oersted field, 
which is calculated from the current-density distribution 
obtained for the nanopillar contact geometry including 
the leads. The Oersted field distributed in the disk does 
not have the circular symmetry. This is due to the asym- 
metric input-output arrangement of the technical current 
(see Fig. [j}. The maximum field value at the perimeter 
of the disk is about 41 mT on one side and 33 mT on the 
opposite side. 

Fig. [3] shows the dependence of the simulated GMR on 
the applied magnetic field for both sweep directions. The 
observed dependence is in a good qualitative agreement 
with the experimentally measured GMR curve for sweep 
2. Since the GMR curves are symmetric in the direction 
of the field sweep we only discuss one sweep direction 
(from positive to negative). For field values larger than 
200 mT a low resistance region is found, in which the 
magnetization in both disks is mostly aligned with the 
external field. At about 190 mT, a vortex nucleates in 
the bottom disk, leading to a jump in the resistance 
which then further increases linearly with the decrease 
of the external field. The vorticity of this vortex is 
determined by the direction of the Oersted field. When 
the external field approaches zero a vortex with the same 
vorticity nucleates in the top disk, which is reflected 
in the low resistance value in this regime. With the 
further decrease of the external field both vortices are 
moved towards the perimeter of the disks. The vortex 
of the upper disk is expelled at about — 320 mT and 
consequently the resistance increases. The resistance 
reaches its minimum above —400 mT when the vortex 
of the bottom disk is expelled as well. The parallel 
orientation of the two magnetizations is restored. The 
"wiggles" in the resistance accompanying the nucleation 
of the first vortex are caused by excess energy. In the 
experiment, this energy is dissipated almost instantly 
and thus they are not observed. 

During the measurement depicted in Fig. [2ja) in ad- 
dition to the dc voltage we measured a spectrum at 
each field value. The star symbols in the graph mark 
the presence of peaks in the spectrum representing ex- 
citations of the magnetization. Obviously, the single 
vortex state is excited by the current. Fig. |2jb) dis- 
plays the spectral power density at each field point. In 
the first field range from 202 to 182 mT the spectra 
are broad (730 — 1740 MHz) and show multiple peaks of 
different amplitudes. As the field is lowered, the fre- 
quencies mostly red shift. At the next field value all 
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secondary peaks have vanished. Starting from 177 mT 
we measure single peak spectra with a peak at initially 
1.275 GHz, which continues the red shift observed before 
until 136 mT, where the frequency reaches its minimum 
of 1.218 GHz with a FWHM of 3.2 MHz. Afterwards 
the peak frequency increases until 106 mT to a value of 
1.356 GHz. 

If we reverse the current sign the resistance profile indi- 
cates that we again have a vortex in a similar field range, 
while the other layer has homogeneous magnetization. 
But the rich dynamics observed with the previous cur- 
rent setting do not occur. In order to rule out thermal 
activation of the vortex motion we prepare the single vor- 
tex state by applying —17 mA and sweeping the applied 
field from positive saturation to 151 mT. After reducing 
the current stepwise to zero while keeping the field con- 
stant we measure a current loop at this field from to 
—20 mA, then to +20 mA and finally back to mA. Only 
for negative current we see clear excitations. Therefore, 
thermal activation cannot explain the excitations, which 
must be caused by spin-transfer torque. 

These results lead us to the following conclusions: 
Since both disks have a geometry, which at low effec- 
tive fields favors the vortex state, and due to the match 
of the resistance profiles obtained by micromagnetic sim- 
ulations to the measured ones we have assured that the 
investigated state comprises one vortex and one homo- 
geneously magnetized disk. The homogeneous disk is 
aligned with the field, therefore its dipolar field reduces 
the applied field to sufficiently small effective values for 
the vortex to occur in the other disk. No significant in- 
plane stray field components are expected to arise from 
the vortex, thus the homogeneous disk is subject to the 
full applied field. If the observed excitations took place 
in the homogeneously magnetized disk, we would expect 
significantly higher frequencies than measured. 

Another very interesting aspect is the V-shaped de- 
pendence of the frequency on the field (Fig. [2jc)). The 
most red shifted peak measured occurs at an applied field 
of 136 mT. This value is close to the micromagnetically 
computed dipolar field that the spacer facing part of the 
vortex disk is exposed to. If we assume that at this point 
the dipolar field and the applied field cancel each other, 
the measured frequency behavior can most easily be ex- 
plained: When the vortex enters the disk the effective 
field is roughly 60 mT. The magnetostatic energy min- 
imum, which is the center of the vortex core motion is 
shifted towards the disk rim. As the field is lowered the 
center of motion moves towards the disk center where 
the frequency reaches its minimum value. By further 
decreasing the applied field the average vortex position 
is shifted towards the opposite rim as the effective field 
tunes down to about — 50 mT. (At lower fields the sam- 
ple enters the double vortex state.) The shifting of the 
average vortex position is connected to a shape-change 
of the potential in which the vortex moves. This change 
of potential shape causes the gyration frequency to first 
decrease on the way to the disk center and afterwards in- 



crease again while the vortex is shifted towards the rim. 

We therefore conclude that the vortex is excited and 
this excitation is caused by spin torque. The vortex 
moves around an average position, which is shifting 
through the disk as the magnetic field is swept. The cur- 
vature of the potential at its minimum is connected to 
the restoring force. With changing field, that curvature 
also changes. This is reflected in the observed V-shaped 
dependence of the frequency on the field. 



IV. VORTEX GYRATION UNDER INFLUENCE 
OF A HOMOGENEOULSY POLARIZED 
CURRENT 

The excitation of the gyrotropic motion of a vortex by 
a spin-polarized current is conventionally described by 
extending the Thiele equation with a spin torque term 
[l3l H4] . However, this approach denies the possibility to 
excite a vortex with a current that is homogeneously in- 
plane polarized [TS] . Khvalkovskiy et al. [TS] pointed out 
that an inhomogeneous polarizer can supply the vortex 
with energy. Here, we will propose another mechanism 
by showing that even a perfectly homogenous polarizer 
can pump energy into the gyrotropic motion of a vortex 
via spin torque, if the angular spin torque efficiency func- 
tion is asymmetric. We start with the Landau-Lifschitz- 
Gilbert equation with spin torque 
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where m and p are vector fields of unit length denoting 
the normalized magnetization and the polarization, re- 
spectively, a is the Gilbert damping parameter, 7 the 
gyromagnetic ratio and j the current density. The effec- 
tive field is defined as 
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where W is the energy density of the disk and M s is 
the saturation magnetization. The spin torque coeffi- 
cients are a = jh/ 1 (2^M s eL) with e being the modulus 
of the electron charge and L the disk thickness. The 
spin-transfer torque efficiency function g((3), where j3 is 
the angle between p and m (cos(/3) = m • p) according 
to Slonczewski [5] reads 
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with P £ [0, 1] being the degree of current polarization. 
It is related to the difference in spin dependent resistances 
R+, Br as P = (R~ - R + )/(R+ + R~). The parame- 
ter A is defined according to A = yj AG{B+ + R~)/2, 
where A is the contact cross section and G is the spacer 
conductance as defined in [8]. Thus A describes the mis- 
match between spacer and ferromagnet resistance and is 
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therefore related to spin accumulation at the nonmagnet- 
ferromagnet interfaces. Let p from now on be homoge- 
neous and constant in time. By combining equations (1) 
and (2) we get the power density absorbed by the disk: 
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The total power is obtained by an integration of equation 
(4) over the disk volume. Let the cylinder axis of the disk 
be the z-axis while the disk is parallel to the x-y-planc. 
The direction of m is then described by the two angles 
<p(x, t) and 6*(x, t) where ip is the azimuthal angle (p = 
is the cc-direction) and 9 is the angle between m and the 



m = sin(#) cos(Lp)e x + sin(6>) sin(p)e y + cos(9)e z . (5) 

We switch to cylindrical disk coordinates (radius p, az- 
imuth x, height z). In order to carry out the integration 
we assume that the magnetization distribution obeys the 
conditions: 
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(6b) 
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a and Xv ar e the polar coordinates of the vortex core, 
respectively. It it assumed that the magnetization dis- 
tribution is independend of the z-coordinate. The condi- 
tions above state that m(a, Xv) is related to m(a, 0) by a 



rotation around the z-axis. They are for example fulfilled 
by the well known double vortex ansatz [16j (using the 
notation of [T7]) 
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where R is the disk radius and 9 is some bell shaped func- 
tion like, for instance according to the ansatz of Thomas 
and Feldtkeller HE] 



= arccos [± exp(— 2« 2 r 2 )] 
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with 



r := v 7 (pcos(x - Xv) - a ) 2 + (/?sin(x - Xv)) 2 - 

In equation (9), the + refers to a positive, the — to a 
negative core polarity and k is determining the core ra- 
dius. In order to obtain an expression for the energy the 
vortex absorbes due to spin torque during one revolution 
when moving on a constant orbit radius a at constant 
frequency u) one has to evaluate the expression 
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Chosing p T = (1,0,0) equation (10) can be transformed 
to 
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with g(A, sin 6* cos ip) — g(A, cos /3) := <?(A,/3). 
k T (xt,) = (sinx«,cosxt),0) is a vector related to the ac- 
tual vortex position. To proceed further it is useful to 
express g in terms of the parameter £ := A 2 — 1. £ = 
now corresponds to the symmetric case where we perform 
a Taylor expansion in £: 

= \ + l ^P l Ptf: (-(co^ + l)*;)" (12) 

This expression only converges for all angles if — 1 < £ < 
1. Keeping the first order (n = 0) term in £ and inserting 



it into equation (11) yields: 

^L-iH 5 ?/^ <13) 

The formula clearly shows that only the core region con- 
tributes to the energy uptake. The sign of the integral 
is independent of core polarity and chirality because vor- 
tices of opposite polarity revolve at different sense of ro- 
tation. The integrand in equation (13) is sensitive to the 
core shape. Since moving vortices are accompanied by a 
dip in the m z -magnetization [19j [20] , we take a magneti- 
zation distribution of a moving vortex obtained from an 
OOMMF-simulation for computation of expression (13). 
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FIG. 4: Time dependence of the total energy of the vor- 
tex derived from OOMMF simulations. A vortex initially 
moving on an intermediate radius is relaxing while spin po- 
larized currents of four kinds are applied. Each simula- 
tion has the same initial state. Thin sinusoidal lines: vor- 
tex energy. Thick lines: Fit of exponential decay function 
y — Aexp(— x/t) + yo. Black: electron flow favors paral- 
lel alignment and A = 1.41. r = 2.4257 ± 0.0008 ns. Red: 
electron flow favors antiparallel alignment and A = 1.41. 
r = 2.0326 ± 0.0005 ns. Blue: electron flow favors parallel 
alignment and A = 1.0. r = 2.2236 ± 0.0005 ns. Green: 
electron flow favors antiparallel alignment and A = 1.0. 
r = 2.2087 ± 0.0005 ns. Inset: All four curves for t < 20ns. 

By comparison to the energy dissipated during one pe- 
riod 
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the critical current for steady gyration is estimated to 
about 3.36 x 10 9 A/cm 2 (for P = 0.85, A = 1.41 and 
L = 30 nm). The vortex absorbes energy for negative j 
(meaning that electrons flow from the polarizer into the 
vortex) and positive £ corresponding to A > 1. If ei- 
ther the sign of j or £ is flipped, the term will lead to 
an additional damping force. The current density es- 
timated above is one order of magnitude higher than 
usually applied currents. This suggests that the effect 
is not the only cause for steady gyration. However, 
simulations clearly confirm it's existence (Fig. [4j. Wc 
perform four simulations (by OOMMF, lateral cell size 
2nm, one z-node, disk diameter 150 nm, disk thickness 
30nm, ^ M S = 2.14T, K 1 =48kJ/m 3 ) with the same 
initial magnetization distribution which is a vortex mov- 
ing on an orbit with an intermediate radius. The ho- 
mogeneous polarizer points into the x-direction. A cur- 
rent of density 1.5 x 10 8 A/cm 2 is employed. In the first 
simulation (black lines) the current sign is chosen such 
that parallel alignment is favored, while A = 1.41 cor- 
responding to £ = 0.988. According to equation (13) 
in first order in £ the spin current must lower the effec- 
tive damping in comparison to simulation three (blue), 



where we have the same current but £ = 0. Simulation 
two (red) is the same as simulation one, but this time the 
electron flow is reversed so that the current now favors 
the antiparallel configuration. In that case we expect the 
strongest damping. Simulation four has the same param- 
eters as three, except that the current sign is negative 
this time. The damping should be the same as in simu- 
lation three. As is shown in Fig. [4] all cases qualitatively 
reproduce the analytical findings: From a fit of an expo- 
nential decay function to the data, damping parameters 
are extracted. Indeed, simulation one exhibits the small- 
est damping corresponding to the largest time constant 
r rj 2.426 ns, simulation two has the largest damping 
(t s» 2.033 ns) and simulation three has an intermediate 
value as expected (r ss 2.224 ns). Simulation four yields 
r w 2.209 ns which is very close to the result of three as 
it should be. In realistic cases A can be much larger. 
For the Fe/Ag/Fe(001) system at low temperatures we 
have an experimental value of A = 3.4 obtained from 
spin-torque measurements leading to £ = 10.6 [S] (and 
references therein). Equation (12) converges uniformly 
for all f3 £ [0,7r] only if |£| < 1, but was well suited 
to calculate the derivative (13). However, it fails if one 
wants to compute the energy contribution for larger |£|. 
In this case one can expand in cos /?: 

PA 2 00 /I - A 2 \ ™ 

n— ^ ^ 

This series uniformly converges on j3 £ [0, tt] for every A. 

V. SUMMARY AND DISCUSSION 

In the first part we presented measurements of excited 
states of a nanopillar containing two ferromagnetic disks. 
One of the disks is in a vortex state while the other is 
homogeneously magnetized. We deduced that the mea- 
sured peaks correspond to spin torque driven vortex mo- 
tion. Several micromagnetic simulations indicate that 
the vortex resides in the thicker bottom disk. We see 
the excitations for electron flow from the bottom to the 
top disk. This would mean that the excitation of the 
vortex is caused by electrons which are reflected by the 
top layer. It is not clear why a direct electron flow from 
the top disk into the vortex cannot excite it. The exper- 
iment in [5] which also has a vortex and a homogeneous 
polarizer shows excitation only for electron flow into the 
vortex. On the other hand, the previously reported ex- 
periment which also involves the Fe/Ag(001) system [TU] 
is in accordance with the presented data. The V-shaped 
frequency vs. field dependence of the vortex excitations 
has been explained by the in-plane motion of the vor- 
tex potential minimum in the disk and the accompany- 
ing shape-change of the potential under the action of the 
field. 

In the second part we have shown that the angular 
asymmetry in the spin-transfer torque efficiency function 
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enables the vortex to gain energy even when the polarizer 
is perfectly homogeneous. The same effect can also 
lead to an additional damping depending on the signs 
of current and asymmetry parameter £. Although the 
OOMMF simulations performed so far are too simple to 
claim capturing the complex behavior of a thick mag- 
netic disk, they clearly show the expected qualitative 
behavior. Nevertheless more advanced simulations will 
have to be performed. This novel mechanism is expected 
to be infuential in the Fe/Ag/Fe(001) system, especially 
in the experiment reported by Lehndorff et al. [10] for 
two reasons: (i) The Fe/Ag/Fe(001) system exhibits a 
rather strong £ = 10.6 and (ii) the polarizing layer is 
laterally extended and thus much less susceptible for 
nonhomogeneous magnetization distributions. 
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